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Abstract 

It has been shown in the previous paper that the metric in the semiclassical 
region of the collapse spacetime is expressed purely kinematically through the Bondi 
charges. Here the Bondi charges are expressed through this metric by calculating 
the vacuum radiation against its background. The result is closed equations for 
the metric and the Bondi charges. Notably, there is a nonvanishing flux of the 
vacuum-induced matter charge. 
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In ref. [1], a new approach is proposed to the problem of backreaction of the Hawking 
radiation. For the spherically symmetric collapse of a compact matter source it has been 
shown that, in the semiclassical region [1] of the expectation-value spacetime, the gravity 
equations close purely kinematically leaving the arbitrariness only in the data functions. 
The data functions are two Bondi charges appearing as coefficients of the expansion of 
the metric at the future null infinity I'^ 
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Here ci is a constant [1], and, for the model of the vacuum considered below, ci = 0. The 
notation in eqs. (1), (2) and below is conventional for spherical symmetry and is the same 
as in ref. [1] with the following exception. The retarded time normalized at the future 
null infinity is in the present notation just and the retarded time u~ counted out by 
an early falling observer [1] will here be denoted as 

u- = U{u) , ^ = U{u) . (3) 

The advanced time v that also figures below remains normalized at the past null infinity 
X~ . The partial derivative du is the derivative at fixed v. 

As shown in ref. [1], the metric in the semiclassical region depends only on the Bondi 
charges }A{u) and Q'^{u). On the other hand, the Bondi charges depend only on the 
metric in the semiclassical region, this dependence being already a subject of the quantum 
dynamics. Having both dependences obtained, one gets self-consistent equations for the 
Bondi charges and, thereby, for the expectation value of the metric in the semiclassical 
region. 

The calculation thus falls into three stages. The first stage: solving the kinemati- 
cal equations for the metric in terms of the Bondi charges is accomplished in ref. [1]. 
The purpose of the present work is the second stage: calculation of the vacuum radia- 
tion against the thus obtained gravitational background. This should express the Bondi 



charges through themselves, and there will remain only the last stage: solving the resul- 
tant self-consistent equations. 

Both Bondi charges A4{u) and Q^{u) appear in the expansion of the flux component 
of the energy- momentum tensor at I'^: 
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Here the integration is over the unit 2-sphere S entering the product X+ = (li-axis) x S. 
Since the collapsing matter source is assumed having a compact spatial support [1], the 
T^^ in eq. (4) is the energy-momentum tensor of the in- vacuum, T^"^. The quantity (4) 
with T^^ = is the only one that needs to be calculated. There is more than one 
way to do this calculation in semiclassical theory. The present calculation uses the WKB 
technique along the lines of refs. [2,3]. 

To model the vacuum, the simplest quantum field will be chosen: a massless scalar 
field satisfying the equation 

(□-ei?)* = 0, e = const. , (5) 
and having the energy-momentum tensor 
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The calculation will be carried out with an arbitrary ^ but the case of interest is ^ = 1/6. 
The latter is important because the locality of the trace of T^c is assumed in ref . [1] . 
The T^^^ at X"*" can be obtained by a direct averaging of the operator (6): 
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and here one of the main points emerges. The subtraction in eq. (7) is necessary because 
there is a noise at X+ even when the background field is absent. The principal requirement 
that the subtraction term must satisfy is its locality in the background field. In the 
effective- action technique, this requirement is secured by the locality of the counterterms. 
In the present technique, it should be secured by the choice of the quantum state in 



the subtraction term. The definition of this state should not involve the quantum-field 
operators in the future or past of the observation point of T^^. With the observation 
point at T"*", the normalization state is the out- vacuum. Below, the notation is introduced 



(X) — (in vac|X|in vac) — (out vac|X|out vac) 



From eq. (6) one obtains 
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Expanding $ in spherical harmonics, one can replace the quantum field in 4 dimensions 
with a sequence of quantum fields in 2 dimensions. The fields I — 0, 1, . . ., are 
defined on the Lorentzian subspace of a spherically symmetric spacetime and, there, satisfy 
the equation 
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with A the D'Alembert operator in this subspace. In terms of one obtains from eq. (9) 
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The expansion of at X+ is readily obtained from eq. (11). One needs the retarded 
solution which is 

l/(/ + l) 
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Here ^i{u) are the data for at X"*". Inserting this expansion in eq. (12) and using eq. 
(4) one obtains finally 
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The flux of in eq. (15) appears as a total derivative (candidate for discarding) but 
appearances may deceive. 

Eqs. (14), (15) leave one with the expectation values of the form 
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where I> — V^du) is some (retarded) linear operator. Let ipl^{x), a — £in, be the solution 
of eq. (11) that asymptotically at becomes the eigenfunction of the energy operator 
with the eigenvalue £in. Let ipi^^{x), A = Sout, be the solution defined by a similar con- 
dition at T+. Normahzed with the aid of the inner product, the ■0;a(^) '^lA'i^) make 
two bases of solutions of eq. (11) (for each Z), related by the Bogolyubov transformation 

^uix) = a^Aa^llix) + Wr(x) . (17) 

(Summation is assumed over the repeating a or/and A but not I.) The dagger designates 
complex conjugation, and "the complex conjugate quantity" will be abbreviated as c.c. 
Expanding the quantum field '^i{x) in any of the bases and calculating the difference in 
eq. (8) one obtains 

{(V^iix))") = -(P^/r(x))(PV^/r(x))/3i„ + c.c. , (18) 

= i/ d'x (')^^/'5(a)(V,a)(C V''V^,T(x) -^rr V^Cl:^)) • (19) 
Here d?x ('^^g^/'^ is the volume element of the Lorentzian subspace, and the equation 

S : a{x) = (20) 

with V(T past-directed describes an arbitrary Cauchy surface E in this subspace. 

To see that the expectation value (18) with x e X+ is causal, i.e., does not contain 
the background field in the future of the observation point, note that it can be expressed 
through the commutator function of the field 

G,{x,x') = \[^i{x),^i{x')] = -iC(a;)C\^') +C.C. = -i^,T(a;)^/r(a;') +C.C. . (21) 

Namely, 



+ (PC V"' ^'t - V^-^f i^[^r')T^Gi{x, x')] . (22) 



With the observation point x at X+ at a given value of u, choose the Cauchy surface S 
as shown in Fig.l, E = E1UE2. The contribution of E2 drops out of (22) because the 
commutator function vanishes when x' is outside the hght cone of x. Then shift Ei to 
T+. The result is 
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Here iIjia{u) and il^la{u) are the data at X+ for the basis functions ■0£4*(^) V'za(^)- 
Only these data and only in the past of the observation point u are needed to calculate 
the expectation value (23). 

The WKB solution for V'zal^) boils down to the solution for the geodesies in the 
background metric. The background metric is the sought for expectation value of the 
metric — the final goal of the calculation — but kinematics gives it in terms of arbitrary 
Bondi charges [1]. One needs to solve for the congruence of the null geodesies that start 
at T~ with one and the same energy and one and the same angular momentum L = hi. 
It suffices to know which of these geodesies come to the point w of X"*", and what then is 
their energy with respect to the KiUing vector at T+. Denote this energy as 

e+ = e+{%£in,L) . (24) 

One can put the question differently. Consider the geodesic that comes to the point u of 
X"*" with the energy and angular momentum L, and trace it back to X^. Then what 
is its energy with respect to the Killing vector at X~? The answer is given by eq. (24) 
solved with respect to Sin- The result is 



. T<, ^^^^ 

£in3 , L<Si^H[U)- . 

du du 
Two functions of the background metric govern this behaviour. One is the U{u) of eq. 

(3), and the other one, H{u), may be interpreted as defining the variable height of the 

centrifugal barrier. Both are expressed through the Bondi charges: 

M{u) + J9M^(u) - 8Q^(u) 
H-\u) = 2 — — ^ — ^ ^ 3 , (26) 

(2>M{u) + ^9M\u) - SQ%u)) 
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k[u) = > u>Uo. (28) 

{M{u) + ^JM^{u) - Q^u)) 

Here uq labels the radial light ray tangent to the apparent horizon [1]. Expression (28) is 

invalid for early u, and the range of u for which eq. (25) holds \su> uq-\- 0{M.) but this 

is the range in which one needs to calculate the radiation. At earlier li, the radiation is 

negligible. The geodesies whose and L do not satisfy either of the inequalities in eq. 

(25) do not come to the point u of X+. 

The derivation of the result above will here be omitted. It will only be mentioned that 
the geodesies corresponding to the second line of eq. (25) turn (in r) three times. The 
respective particles start at X~ as incoming and turn the first time before the black hole 
has formed. They are outgoing already when the collapsing mass comes and turns them 
back. They cross the apparent horizon and continue falling down but, before they reach 
small r and even before they cross a possible event horizon [1], gravity weakens and lets 
them go. They cross the apparent horizon the second time, next turn the third time, and 
go out to . The third turn is a quantum effect. 

The congruence of geodesies considered is hypersurface orthogonal, and the function 
(24) defines the phase of '^/^(a;) at X+. Hence one obtains the data at X+ for the basis 
functions: 

C(ii) = \e{e,^HU - I) exp(-i£inC/) + e{l - Si^H) exp(-i£inK)l , (29) 

^,T(m) = ^^J_ exp(-i£out^). (30) 

The expectation values (23) will thus be expressed entirely through the functions H{u) 
and U{u). It is important that, even kincmatically, the Bondi charges are not completely 
arbitrary [1]. By their properties, H{u) and n{u) in eq. (27) are macroscopic quantities 
whereas the derivatives of these functions in u are microscopic quantities. In the notation 
ofref. [1], 

\0\<H,^-<^ ^H^O, ^--O, (31) 
K, \U\ du dun 



7 



il = , u>uo (32) 

where O vanishes when the quantum parameter tends to zero. 

In the upper hmit of the integral in eq. (23) one may exploit condition (32): 

U(u + 0) <U{u) (33) 

and set U{u+0) — 0. In view of eqs. (27) and (31), this is equivalent to setting u+0 — oo, 
and then the integration by parts proves that the two terms of the integral are equal. In 
the product of two functions (29) that figures in eq. (23) only one diagonal term survives 
and takes the form 

^ll^'i^) = / 4^ ^(^in^^ - ew{ieUU - U')) (34) 

™ 

because, for u > u', HU < H'U' . The contribution of the other diagonal term to the 
expectation values vanishes because the contribution of the term (34) vanishes when 
K = 0, and expression (23) vanishes when = ■0°"*. The cross terms vanish because 
they imply U > \0\ or U' > \0\. 

It follows from eqs. (23) and (34) that 

oo °° 

Y.p,{l){{v^i{u)f)=p,m{v^i{u)f) + / dlp^{l){{v^i{u)f) (35) 

/=0 '=0 

where Pi{l) is any given polynomial, and P2{1) is some other polynomial. Doing the sum 
over / first and the sum over e^yx last, one obtains each contribution to (35) as a spectral 
integral over the energy £out, the spectral function being a combination of the functions 

oo oo 

In{£out,u) = j dSinU (iSinU) j (iw' £in^/' exp(i(f] - f]')) , (36) 
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n -n' = £j„([/ -[/') + £„^t(M - k') ■ (37) 

Specifically, 
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In eq. (36) introduce the new integration variables 



y = 6inU- , X = £inf/'^ (40) 



to obtain 



oo 



/n(£out, u) = J dy {iy f J exp(i(l] - ^l')) , (41) 



Next use eq. (27) to write 

u 

if = exp(- J duKj , (43) 

— oo 

u' u" 

U'-U^ J du" exp(- j du r) (44) 

u — oo 

and integrate in u" by parts: 

y du K du \K, du K J J n 

-d — - —(- — -\ \-tf' (45) 
\ ^ du'K'^ du'\K'du'K')^"') k' ' ^ ^ 

Owing to condition (31), all the corrections with the derivatives of k are negligible both 
in (42) and (45): 

w^l, £in(C/ -U')^x-y. (46) 

There remains to be expressed through x and y the difference {u — u') in eq. (37). For 
that one has the equation 

ln^= 1 du"K"{l-'^-\ (47) 

X J \ du" k" ^ ' 

in which the last term is negligible. This equation can be solved by expanding in the 
derivatives of k: 

.(«'-«)= 1. ^ + i f ^i) + i f f f i^i) ) 1„» 1^ + . . . (48) 

X I \du Kj X 6 \du \Kdu Kj J x 

but here the corrections with the derivatives are not unconditionally negligible as they 
are in eqs. (42) and (45). The point here is that the integrals in x and y are cut off by 
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oscillations at both the upper and lower limits so that the integration regions for In x and 
Iny are effectively compact. Then ln{y/x) is bounded, and, in eq. (48), all the terms 
with the derivatives of k, are negligible. As a consequence, the spectral function (41) is 
calculable. Another important consequence is 



(49) 



It was tacitly assumed in the derivations above that conditions (31)-(32) valid at the 
observation point u are valid also at the integration point u'. Eq. (49) proves that this is 
the case. 

Now one comes to the central point of the present consideration. The argument above 
about the effective compactness of the integration regions for In x and In y needs a reserve. 
At the upper limits, the integrals in x and y are cut off by oscillations always whereas at 
the lower limits only when Eout 0- Therefore, the argument may break down for the 
low-energy part of the spectrum: Eout 0. This can be checked. Introduce the notation 

^out 



(50) 



If the corrections in eq. (48) are small indeed, their contributions to the spectral function 
can be calculated by expanding the exponential 



Here P is a series of the form 



i2;,ln — ) explizln — ). 



(51) 



Pfi^,ln-) =yC)(i^)'=fln-)' , k>l, s>k+l. (52) 

V yJ ^ \ yJ 

Insertion of the expansion (51) in eq. (41) yields the following result: 



4(eout,M) = /^"^ («) 



Fn{z) , 



(53) 

F^{z) = e-^'Tin + 1 - iz)r(l + iz) (54) 

where F is the Euler's function. When n > 0, the function Fn{z) and all its derivatives 
are bounded including at 2; — > 0. Then the entire contribution of P in eq. (53) is O and 
is negligible. However, when n < 0, the function Fn{z) behaves as I/2; at 2; — > 0. Then, 
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because, in the series P, the power of (d/diz) exceeds the power of (i^;) at least by one, 
the corrections due to P are even more singular, and, therefore, at £out ^ they are not 
smalL 

Thus, for n > the spectral function (36) is successfully calculated. One has 

n>0: In(eont,u) ^ K''+\u)e~'''r(n + l-iz)r{l + iz) (55) 
with z in eq. (50). Hence one obtains, in particular, the expectation value (38) 



4 7, 27r£out ..^x 
deout -^r; r (56) 
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but does not obtain (39). 

For In with n < one has a difficulty. In order that the function k,{u) could be 
regarded as slowly varying, the operator V acting on ■ip'^"\u) in eq. (23) should be 9„ to 
the power 1 or higher. This suggests the way of overcoming the difficulty. For T> = 
and V — du~^ write in eq. (23) 

u 

V^t^^(M) = J dud^i)^''' , (57) 

— oo 

u u 

j dui)'^'''{u) = j du{u-u)d^^^''' . (58) 

— oo — oo 

Alternatively, write in eq. (36) 

It 

exp(i£inC/) = J a(?i (^isi„f7^ exp(i£inC/) (59) 



to obtain 



u 

/_i(£out,M) = y" du Io{e out, u) exp{i£o^t{u - u)) . (60) 



This calculates /_i through Iq, and for Iq one has the result (55) but the new obstacle is 
that the integral (60) involves u down to u = — oo whereas the result (55) is valid only 
for u > Uq. Indeed, only at late u is k{u) slowly varying by virtue of condition (31). The 
obstacle is not big, however. Expression (55) is inaccurate at early time but, since the 
radiation at early time is negligible altogether, this inaccuracy is unessential. The specific 
form of k{u) at early time is also unessential. It is only important that k{u) falls off at 
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u — oo so that the integral in eq. (43) converges. For the calculational purposes, one 
may just set k{u) = for m < uq, this being equivalent to switching off the background 
curvature at early time. 

Using eqs. (60) and (55), one obtains the expectation value (39): 



de. 



out 



du 



2n sin {sontiu - u)) 



- 1 
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where z — £out//^, and R — k,{u). The integral over ^out can be done: 

4 



(61) 



and finally one obtains 
4 
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du K 
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(62) 
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du Kj—Tilnl / duK,]+0{l) 



u 

j duRj 



oo . 
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This is the message of the present paper. The operator averaged over the in 



x+ 



obtained as 



vacuum or out-vacuum is infrared-divergent. The expectation value ( 
the difference in eq. (8) is finite but growing at late time. Its first derivative contributes 
significantly to the fiuxes (14)-(15): 



=0 (47r)2 2 

whereas the second derivative is already negligible: 



^ ^k{u){1 + 0) 



(64) 



(65) 



=0 "-dw 

This mechanism of emergence of the vacuum fluxes is familiar. In the effective-action 
technique they all emerge as total derivatives of growing vertex functions [4] . 

An alternative way of calculating expression (61) is introducing the integration vari- 
ables 

7 — £out('^ ~ u) , a — k{u — u) , u > Uq (66) 

and doing the integral over a first. It works also when calculating the / > contributions 
to the expectation values. Using eqs. (57), (58) one obtains 
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X < 



-(-l)tln( J duR) +0(1) , 

— oo 

I 0(1) , 



n even 



n odd , 



(67) 



Jdir{i J du^i{u) 



1 n\ 
(47r)2 n + 1 



. n 

-1)^ 



i/"+^(ii)K"-^(ii) 



n even 



X < 



(-l)^27r(y dunj +0(\n(^j duRj) , n odd . 



(68) 



Here the powers of growth are different for even and odd powers of / but in all cases this 
growth is insufficient. Since, in eq. (15), on the terms (67) there acts du and on the 
terms (68) the contributions of all of these terms are negligible. Only the s-mode 

contributes to the flux of the charge Q^. 

The I > modes contribute only to the flux of the gravitational charge M. via 

2 1 



j dir {{d^^iiu))' 



(47r)2 n+l 



X j d£out ^out [-(-i)"e-''T(n + 1 - iz)r(l + i^) + c.c. 





(69) 



The WKB technique is inaccurate for I of order 1 but the contribution of the s-mode 
is unambiguous and so is the contribution of the / ^ 1 modes. The latter is given by 
the highest power of I in the polynomial P2(0 ®<1- (35)- Retaining only these two 
contributions one obtains 



3 

out 



(70) 



The contribution of the s-mode describes correctly the low-energy behaviour of the spectral 
function, and the contribution of the I » 1 modes describes correctly the high-energy 
behaviour. The inaccuracy at intermediate energies is a question of the grey-body factor. 
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The final result is the following set of equations for the Bondi charges: 

"^"•^ = "'^""^ {^^Yo mu)<^)f) , (71) 

9uQ' = ^«:H(l-40. (72) 

The contribution of the / ^ 1 modes to the total energy flux (71) is by an order of 
magnitude less than the contribution of the s-mode because the high-energy part of the 
Planckian spectrum is exponentially suppressed. The flux of the charge depends on 
the value of ^, and there are values for which it is zero or negative but, if ^ 7^ 1/6, the 
present calculation is inconsistent because the background metric of ref. [1] is invalid. 
The nonvanishing flux of the charge is a surprise. 

In conclusion, a failure of the 2-dimensional effective action is worth mentioning. At 
^ = 0, to the s-mode *o there corresponds the effective action in the Lorentzian subspace: 

(^^^vac = / d'x (V (73) 

so that 

(%v^ (v^;:v^ ^ ^ ^ " °^ ^ ^ ^ 

(with the retarded current on the left-hand side). Eq. (74) is to be compared with eq. 
(12) at ^ = and / = 0. Only to the leading order in 1/r do these expressions coincide. 
The effective action (73) should, therefore, give the correct energy flux for the s-mode, 
and it does. But the flux of the charge is not contained in this action even for the 
s-mode and even at ^ = 0. This may explain why the 2-dimensional models of the effective 
equations miss the backreaction of radiation. The 4-dimensional effective action should, 
of course, reproduce the present results in full but here it will not be considered. 

The equations above for the Bondi charges close. Thereby, the expectation-value 
equations for the metric close already at the level of functions of one variable [1]. The 
solution will be reported. 
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Figure caption 



Fig.l. Choice of the Cauchy surface E for the inner product, E = Ei [j^2- The division 

into El and E2 is specified by the location of the observation point u at X"*". The 
event horizon EH if any is in the future of E. 
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